Iterated fully coordinated percolation on a square lattice 
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We study, on a square lattice, an extension to fully coordinated percolation which we call iterated 
fully coordinated percolation. In fully coordinated percolation, sites become occupied if all four 
of its nearest neighbors are also occupied. Repeating this site selection process again yields the 
iterated fully coordinated percolation model. Our results show a substantial enhancement in the 
size of highly connected regions after each iteration (from ordinary to fully coordinated and then 
to iterated fully coordinated percolation). However, using Monte Carlo methods to determine the 
static critical exponents and normal mode analyses to determine the dynamic critical exponents, 
our results indicate that the ordinary, fully coordinated, and iterated fully coordinated percolation 
models belong to the same universality class. This is in contrast to our previous results wherein 
a different universality class dynamically (but not statically) for fully coordinated percolation was 
indicated, though the differences in the dynamic exponents were small. A possible cause might be 
the different methods of generating clusters, i.e., in this work clusters were generated statically in 
square lattices while in our previous work clusters were grown dynamically, for studying the dynamic 
critical exponents. 

PACS number(s); 64.60. Ak, 64.60. Fr, 05.40.Fb, 05.70.Fh 



I. INTRODUCTION 

Iterated fully coordinated percolation (referred to as 
IFC percolation from this point on) is an extension to the 
ordinary percolation and the fully coordinated percolation 
(referred to as FC percolation hereon) models. These 
are well-known problems of geometrical phase transition 
where global connectivity is singular as a function of lo- 
cal connectivity (see, for a review, Stauffer and Aharony 
[^). In ordinary site percolation on a square lattice, sites 
are randomly occupied with an independent occupation 
probability p. Two nearest neighbor sites are considered 
connected if they are both occupied. A collection of con- 
nected sites, i.e., every site in the group is connected to 
every other site either directly or through the paths of 
nearest neighbors, is termed a cluster. When p is in- 
creased until it reaches a critical value Pc, an infinitely 
spanning cluster will form for the first time. The square 
lattice shown in Fig. |l](a) is an example of ordinary site 
percolation on a square lattice. Occupied sites are shown 
as circles while unoccupied sites are not shown. 

To perform IFC percolation, first perform ordinary 
percolation. After this is done, determine which of the 
occupied sites have all four of their nearest neighbors 
also occupied. Occupied sites having this property are 
called fully coordinated sites and will remain occupied 
while the other sites will be removed. This process of 
winnowing occupied sites was described in ^ and we 
called the resulting model the fully coordinated percola- 
tion model (FC percolation). The square lattice shown in 
Fig. 0(b) is an example of fully coordinated percolation 
constructed from the ordinary percolation sites shown in 
Fig. 0(a). The gray-filled circles are fully coordinated 
sites. The final step in the process is to again choose 
fully coordinated sites but this time from the resulting 



system constructed from FC percolation. As in the FC 
percolation case, we also remove the sites that are not 
fully coordinated. The solid circles in the square lattice 
shown in Fig. |^(c) are examples of IFC percolation sites 
constructed from Figs. 0(a) and |l|(" ' 
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FIG. 1. An illustration of the process to simulate IFC per- 
colation, (a) First we simulate ordinary percolation on a 
square lattice. Occupied sites are shown as circles while unoc- 
cupied sites are not shown, (b) From the resulting sites from 
ordinary percolation, choose the fully coordinated sites, i.e., 
occupied sites whose four neighbors are also occupied. These 
sites are shown as gray-filled circles, (c) Repeat the culling 
of non- fully coordinated sites one more time. The remaining 
sites are then considered as sites in IFC percolation. These 
sites are shown as solid circles in the figure. 

Notice that this problem is different from bootstrap 
percolation (see, for a review, Adler ||^) where sites of 
less connectivity are iteratively removed until no more 
sites have insufficient connectivities. This is in contrast 
to our problem where we only perform the culling of non- 
fuUy coordinated sites twice. 

In our earlier work |l[], we found that the sizes of 
the fully-connected subregions (with no internal holes) in 
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an incipient infinite cluster in FC percolation are much 
larger than those found in ordinary percolation. In this 
study, we find this same increase in the sizes of the fully- 
connected subregions to continue when the full coordina- 
tion requirement is iterated one more time. The regions 
of high local connectivity are generally interesting be- 
cause they would seriously affect the vibrational modes 
in an equivalent disordered network of springs. In par- 
ticular, it was found previously that they enter into the 
scaling of the normal mode spectrum Q . 

In the same previous work , we also concluded that 
the dynamic critical exponents of FC percolation were 
close to but different from those of ordinary percolation, 
drawing this conclusion from numerical results on the 
FC percolation clusters which were grown from a seed 
site until prescribed sizes were reached. In this work, 
on the other hand, we obtain numerical results on the 
corresponding dynamical critical exponents of FC and 
IFC percolation from the clusters generated statically on 
square grids of prescribed sizes. Surprisingly, our new 
results indicate that the numerical values of these expo- 
nents are very similar for ordinary, FC, and IFC perco- 
lation, and in fact most likely to be the same. Since the 
actual numerical estimates for FC cases are compatible 
within the error estimates of both works, this might sim- 
ply mean that we drew an overzealous conclusion in our 
previous work. However, we cannot rule out the possi- 
bility that the ensembles of clusters used in and the 
current work may contain subtle differences which in fact 
alter dynamic universality classes. 

The rest of this paper is divided into four sections. In 
Sec. |l| we present numerical results and calculations of 
the static and dynamic critical exponents of IFC perco- 
lation. The geometry of clusters is studied in Sec. Ill 
with a particular focus on the sizes and distributions of 
the highly connected subregions. The difference between 
the methods of generating the cluster ensemble in [Q and 

The final section 




the current work is discussed in Sec. ^ 
gives a summary of our results. 



FIG. 2. Susceptibilities x against probability p in IFC per- 
colation on square lattices of sizes are shown. The squares, 
circles, diamonds, and triangles are for L — 256, 512, 1024, 
and 2048, respectively. Each data point in the plot is an aver- 
age over 1000 realizations. Shown as the inset is the finite-size 
scaling of the susceptibility peaks by plotting the four peaks 
against L on a log-log scale. Note that L is a dimension- 
less quantity corresponding to the length, in multiples of the 
lattice constant, of the side of the square lattice. 



The critical occupation probability, Pc, is determined 
from the plot of susceptibility, x, against p. This plot 
is shown in Fig. |[ Notice that because of finite size ef- 
fects, instead of a diverging x at Pc, we see clearly defined 
peaks at p values near the true pc- The lattice sizes, L^, 
we used in our simulations were L = 256, 512, 1024, and 
2048, and each data point in Fig. || is an average over 
1000 realizations. According to the plot, the suscepti- 
bility peaks are pegged at Ppeak = 0.946 for the larger 
lattices. Since the larger the lattices the more our simu- 
lations mimic the asymptotic limit, we can conservatively 
assign the value for the critical occupation probability to 
be atpc = 0.946 ±0.001. 



II. STATIC AND DYNAMIC CRITICAL 
BEHAVIOR 



We have determined the static and dynamic critical 
exponents of IFC percolation by simulations on a square 
grid. The culling of non-fuUy coordinated sites occurs in 
two stages. After sites in the grid are randomly occupied 
with an independent probability p, all the fully coordi- 
nated sites are marked. After this marking process is 
done, all the occupied but unmarked sites are removed. 
From the resulting system of remaining sites, we again 
determine which sites are fully coordinated and remove 
those which are not. Figs. |l| depict this process. To de- 
termine the connectivities among the sites, we have used 
the standard Hoshen-Kopelman algorithm . 



TABLE I. The ratio between static critical exponents 7 
and v, and the fractal dimension df found for ordinary, FC 
and IFC percolation models are shown. The values at the 
bottom row of the table are previously determined values for 
ordinary percolation. The y/f ratio for ordinary percolation 
is not determined in this work. 



percolation type 


7/1^ 


df 


ordinary 




1.888 ± 0.004 


FC 


1.791 ±0.006 


1.885 ±0.014 


IFC 


1.805 ±0.018 


1.857 ± 0.007 


ordinary 
(prior results) 


If 1.7917) " 


§ 1.8958) " 



^see, e.g., Stauffer and Aharony 
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The ratio between critical exponents 7 and v can be 
determined from the scahng of susceptibiUty peak values 
against the corresponding lattice sizes, i.e., x{Ppeak, L) ~ 
jjy/>^_ Using this relationship we find 'j/v for IFC perco- 
lation to be 1.805 ± 0.018. This value agrees nicely with 
the ordinary percolation value which is known exactly to 
be 43/24 « 1.7917. The ratio ^/v can also be used to 
indirectly determine the fractal dimension, df, from the 
known scaling relationship, df — l + ^/2v. The result for 
IFC percolation is df = 1.9G3±0.GG9 while the known ex- 
act value for ordinary percolation is 91/48 « 1.8958. The 
critical exponent 7 can also be determined independently 
from the scaling of susceptibility against the occupation 
probability p, which goes as x{Pj ^ — ^) ^ |p ~ Pc| 
where the L = 00 indicates the scaling to be ideally 
satisfied in the asymptotic limit oi L —^ 00. As in the 
case of PC percolation, the susceptibility peaks are very 
close to p = 1.0. This provides data to the right of the 
peaks to reside in a small interval. We therefore only use 
data to the left of the peaks to determine the scaling im- 
plied. The value we found for 7 in IFC percolation using 
Pc = 0.946 is 7 = 2.398 ± 0.039. This value also agrees 
very well with the value in ordinary percolation, known 
to be exactly 43/18 2.3889. 

The fractal dimension, df, of IFC percolation can 
also be determined directly. At Pc, the size of the 
largest cluster, S, scales with L as S* ~ L'^'f , where 
df = 1.857 ± 0.007. This value is close to the known ex- 
act value for ordinary percolation as well as to the value 
obtainable from our measured 7/j^. 

Listed in Table || are the values for the ratio 7/t/ and 
for df we found from our simulations of ordinary, FC, 
and IFC percolation. As a reference, the previously de- 
termined values of the corresponding critical exponents 
for ordinary percolation are listed at the bottom row of 
the table. 

The values for the spectral dimension, ds, and the walk 
dimension, dyj, can be determined from the dynamic crit- 
ical behavior of IFC percolation. We first construct the 
transition probability matrix W. The elements, Wij, of 
this matrix are the hopping probabilities per step of a 
random walker to hop from site j to site i. The diagonal 
elements Wu are the probabilities for the walker to stay 
at site i at a time step. In the myopic ant rule, all of 
the Wu are zero. In contrast, in the blind ant rule, the 
Wu are nonzero, their values dependent on how many 
occupied neighboring sites the walker can hop into (see, 
for a review, We are using the blind ant rule for this 
work. 

Once the W matrix is determined, the value of dw can 
be estimated directly from its largest non-trivial eigen- 
value, Ai. This eigenvalue satisfies the finite size scaling 
law I 

\lnXi\ w 1-Ai - 5-2/'^^ (1) 

Making use of the Alexander-Orbach scaling law 
ds — 2df/dw, and S ^ L'^f , we obtain 



1-Ai - L-"^-. (2) 

Eq. ^ gives a rather direct way of estimating dw from 
finite size scaling, in a manner that is considerably sim- 
pler compared to methods wherein a detailed knowledge 
of the eigenvectors is required (see, e.g., JTot). 




lattice length, L 



FIG. 3. Shown is the scaling of the largest nontrivial eigen- 
value Ai by plotting 1 — Ai against the lattice length L on a 
log-log scale. The circles, squares and triangles are for ordi- 
nary, FC, and IFC percolation, respectively. Each data point 
in the plot is an average over 1000 random realizations. 

Shown in Fig. ^ is our data plotted in a way that is 
consistent with Eq. (H). For each percolation type, we 
used 13 different edge lengths of the square lattices that 
range from L = 94 to L = 198. Furthermore, for each 
data point, about 1000 random realizations of the W ma- 
trix had been generated by Monte Carlo simulations. To 
numerically determine Ai from each W matrix, we used 
the Arnoldi-Saad method [ pl| . 



TABLE II. Shown are the walk and spectral dimensions for 
ordinary, FC, and IFC percolation. The walk dimensions were 
determined by taking the largest non-trivial eigenvalue of W 
while the spectral dimensions are derived from ds — 2d//du,, 
with df values taken from Table I. The last row in the table 
contains previously known values for ordinary percolation. 



percolation type 


dw 


ds 


ordinary 


2.841 ±0.015 


1.329 ±0.010 


FC 


2.843 ±0.015 


1.326 ±0.017 


IFC 


2.857 ±0.015 


1.300 ±0.012 


ordinary 
(prior results) 


2.87 ± 0.02 


1.30 ±0.02 ^ 



^see Majid, et. al. Q. 
''see, e.g., Nakanishi 



3 



Listed in the second column of Table |l| are the d^, val- 
ues we found following Eq. (^) for ordinary, FC and IFC 
percolation. Notice that there is an excellent agreement 
among the values for all three cases and that they agree 
with the previously known value for ordinary percolation. 

In the third column of Table || are listed the dg values 
calculated using the Alexander-Orbach scaling law. Al- 
though all three estimates from this work agree with the 
previously known value for ordinary percolation within 
the aggregate error estimates of our work and those of 
the previous work, there is a slight internal discrepancy 
between our own estimate of ds for ordinary percolation 
and that for IFC percolation. Since ds is derived from df 
and dw, the slight discrepancy may be attributed to the 
df value we found for IFC percolation by direct simula- 
tion (refer to the third column of Table ^ . We suspect an 
unaccounted systematic error occuring when we perform 
each iteration process. In particular, as we will discuss 
in more details later, the length scales associated with 
fully-coordinated subregions appear to increase at each 
iteration. This suggests that we may have to increase 
the overall lattice length scale correspondingly in order 
to maintain a good estimate of the fractal dimension by 
direct simulation. Since this was not done in this work, 
a systematic error may have crept in this way. 



III. GEOMETRY OF CLUSTERS 

Previously it was found that the number of inte- 
rior sites in the largest cluster in FC percolation is much 
greater than that in ordinary percolation at their respec- 
tive thresholds. In this work we find that this trend 
continues further for IFC percolation. As an illustra- 
tion, consider Figs. ||(a) and ^(b). Shown in Fig. ^(a) 
is the largest cluster from a simulation of ordinary per- 
colation on a 200x200 square lattice at its percolation 
threshold. The darker dotted sites are fully- coordinated 
(FC) sites while the lighter gray dotted sites are partially- 
coordinated (PC) sites. FC sites are those all of whose 
neighbors are also occupied (thus these are similar to, 
though not identical with, the interior sites considered 
in Ql), while PC sites have at least one nearest neighbor 
not occupied. The same shading tone is used in Fig. ^(b) 
except that what is shown is the largest cluster from 
a simulation of IFC percolation at its own percolation 
threshold. Notice the difference in the number of FC 
sites in these figures. Partially coordinated sites domi- 
nate in ordinary percolation while FC sites dominate in 
IFC percolation. 

The average number of FC and PC sites in ordinary, 
FC, and IFC percolation are plotted in Fig. ^. Linear 
least squares fits indicate that in ordinary percolation 
the PC sites are more than 6 times as many as the FC 
sites, while the ratio in FC percolation is about 1.5, and 
in IFC percolation it is only about 0.86. A trend that 
can easily be seen here is that the number of FC sites in- 



creases as we progress from ordinary to FC and then to 
IFC percolation. Although only a local full-coordination 
correlation is imposed, it clearly results in a change on 
the global scale. This may be reflected in some qualita- 
tive differences in a physical property if it is sensitive to 
the size of connected clusters of FC sites. For example, 
if a percolation model is used to describe a disordered 
substrate on which diffusive physics takes place (such as 
in catalysis), what is expected within a given time scale 
(or frequency range) would be directly affected by the 
size of the high-connectivity or FC regions of the sub- 
strate where fast diffusion occurs compared to less well 
connected areas. In such cases, it is expected that there 
would be large differences depending on whether the ap- 
propriate percolation model to use is ordinary, FC, or 
IFC percolation. The latter cases could be a better rep- 
resentation if the culling of partially coordinated sites 
naturally occurs in the particular problem, say, due to 
desorption. 
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FIG. 4. Sample largest clusters resulting from simulating 
ordinary and IFC percolation at their respective percolation 
thresholds in a 200x200 square lattice with free boundaries 
are shown. The darker dotted sites are FC (fully coordinated) 
sites while the lighter gray dotted sites are PC (partially co- 
ordinated) sites, (a) For ordinary percolation, PC sites dom- 
inate. The occasional FC sites are clumped together and are 
somewhat distributed evenly in the cluster, (b) For IFC per- 
colation, the number of FC sites increases. 
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lattice length, L 



FIG. 5. A log-log plot of the number of FC and PC sites 
in the largest cluster at the percolation threshold against the 
lattice length L is shown. The circles are for PC sites while 
the solid circles are for FC sites, both in ordinary percolation. 
The data points for FC percolation are the squares (which are 
PC sites) and the solid squares (which are FC sites). For IFC 
percolation, the triangles are for PC sites while the solid tri- 
angles are for FC sites. Each data point is an average over 
2000 random realizations. 




cluster size, s 

FIG. 6. A distribution of cluster sizes is shown in a log-log 
plot. The circles, squares and triangles are for ordinary, FC, 
and IFC percolation, respectively. Each data point is an av- 
erage over 2000 random realizations. The line drawn is not 
the best fitting line but instead one that follows Nn ~ s~'^ . 
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FIG. 7. Shown are semi-log plots for the number of clusters 
formed by FC sites against the corresponding cluster sizes. 
Each line drawn in each plot does not correspond to the best 
fit line for all available points but only for the four middle 
points. This would give a rough estimate of the average FC 
cluster sizes, (a) In ordinary percolation, on average there 
are 2 FC sites in each FC cluster, (b) In FC percolation, on 
average FC clusters have 18 FC sites, (c) In IFC percolation, 
there are around 72 FC sites in an FC cluster. 

The distribution of cluster sizes (of all sites, i.e., both 
FC and PC sites) is shown in Fig. ^. Because of the fi- 
nite size effects, only the intermediate size ranges (about 
10^ to 10'' sites) are consistent with the cluster size scal- 
ing Ns ^ s~'^ with the ordinary percolation value of 
T = 187/91 « 2.055, where s is the cluster size and Ng is 
the number of clusters having size s. Very small s values 
are clearly out of the asymptotic scaling regime while the 
region of s ~ 0{L'^) is also limited by the finite size of the 
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grid. Nonetheless, it is important to note that there is no 
significant difference between ordinary, FC, and IFC per- 
colation for all s above about 10^. This suggests strongly 
that the asymptotic cluster size scaling (which is known 
to occur for ordinary percolation) also occurs for both 
FC and IFC percolation. 

In making the plot for Fig. ^, we didn't distinguish be- 
tween PC and FC sites in the clusters. If we instead just 
consider FC sites and make a semi-log plot of the number 
of clusters formed by these sites against the correspond- 
ing cluster sizes, we should be able to roughly estimate 
the average cluster sizes formed by FC sites alone. Shown 
in Fig. ^(a) is Ns against s of clusters formed by FC sites 
in ordinary percolation. The line shown is not a best 
fitting line for all of the available points but instead one 
that best fits only the four middle points and is shown 
merely to obtain a crude idea of the typical size scale 
of FC clusters. For ordinary percolation, we found the 
average FC cluster size to be around sq ~ 2 FC sites. 
Shown in Figs. 0(b) and 0(c) are data for FC and IFC 
percolation, respectively. We found that around si « 18 
FC sites comprise the average cluster in FC percolation 
while around S2 ~ 72 FC sites make up the average FC 
cluster in IFC percolation. This means that in order to 
map the FC and IFC percolation problems to the same 
length scale as ordinary percolation, a length rescaling 
by \/si/sq ~ 3 and \/ S2/ sq w 6 must be performed, 
respectively. Thus, although both static and dynamic 
critical exponents seem to remain the same as the full- 
coordination rules are iteratively applied, the local con- 
nectivity greatly changes or, equivalently, the length scale 
which would preserve the local connectivity changes. 

IV. COMPARISON WITH PREVIOUS WORK ON 
FC PERCOLATION 

In this work, we follow the prescription described in 
the Introduction closely in order to generate the ensem- 
ble of FC and then IFC clusters. The only limitation in 
the Monte Carlo portion of our work are that the grid 
we use to generate the clusters must be finite and that 
a particular boundary condition must be selected. Our 
largest grid is 2048x2048 and clearly, this is not very large 
by today's standards. Since we use finite size scaling in 
the determination of both the static and dynamic criti- 
cal exponents, the finiteness of the grid (and cluster) is 
in fact necessary; yet, there is always a possibility that 
an asymptotic regime of large sizes may not have been 
reached. Also, we use the free boundary condition with 
no direct connectivity across the opposite edges of the 
grid, which may accentuate the finite size limitations. 
This type of ensemble is the most straightforward exten- 
sion of the standard Monte Carlo simulation of ordinary 
percolation to the FC and IFC percolation models, and 
may be described as static since there is no dependence 
on the order in which sites are generated or the connec- 
tivities checked. Indeed, the static exponents 7, ^/v, 



and the like, agree closely with the corresponding ordi- 
nary static percolation values, and, if the length scale is 
identified with the grid edge length, then the cluster's 
fractal dimensions should also agree closely with that of 
ordinary percolation. 




FIG. 8. Figure to illustrate cluster geometry difFerences be- 
tween growing a cluster from a seed and generating it in a 
square lattice (static ensemble). The solid circle is the seed 
and the cluster grown from it includes all the circles but ex- 
cluding the gray- filled ones. The growth process stops when a 
specified cluster size is reached. In this figure, the gray-filled 
circles are excluded because the specified cluster size has al- 
ready been reached; however, when generating the cluster in 
a square lattice of the same length scale, the gray-filled circles 
will be included in the cluster. 

In our previous work |Q| , we studied the static and dy- 
namic critical behavior of FC percolation using the en- 
semble of clusters of specified number of sites, say, S' . 
Those clusters were grown from a seed site by applying 
a breadth- first algorithm twice; i.e., first generating an 
ordinary percolation cluster in nearest neighbor shells to 
a sufficiently large size (say, 5S"), and then identifying 
and joining only the fully coordinated sites in a second 
searching process from the same seed and then stopping 
the process when the size of the resulting connected clus- 
ter reaches the predetermined value S' . Now, if we make 
the intermediate cluster of the first search process very 
large and if we do not terminate the second search pro- 
cess arbitrarily when S' is reached, but rather continue 
the search until the natural boundary of the FC cluster 
is found, then the resulting cluster ensemble should be 
exactly the same as that which would result in the static 
ensemble of FC percolation clusters generated on an in- 
finite grid. However, since the first stage is stopped at at 
most 55" and since the second stage is forcibly stopped 
exactly at S", which is a predefined value, we introduced 
a kinetic nature to the cluster ensemble. This may be 
somewhat similar to the differences in the degree of in- 
trinsic anisotropy in static and kinetic ensemble of per- 
colation clusters studied previously (l^Jl^ (called equlib- 
rium versus growing clusters there) . 
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An illustration of the cluster geometry differences aris- 
ing between a static cluster and a growing one is shown in 
Fig. ^. In the static ensemble, every site of the spanning, 
connected component of the ordinary, FC, or IFC per- 
colation models' realizations within a given length scale 
L (defined to be the grid edge length) is a part of the 
cluster. In the growing cluster, on the other hand, only 
those sites which are within a certain path length along 
the nearest neighbor connections from the seed site are 
counted, independent of how close they are by straight 
line distance to the seed. Also, since cluster size S' is the 
control parameter for this ensemble, the length scale L 
is not precisely fixed but rather obtained as an average 
for a given value of S' . Exactly which sites are counted 
in a cluster of length scale L and which ones excluded 
do clearly depend on the particular search order that is 
employed. Thus, typically, for the same L, the number 
of sites of the percolating cluster is smaller in the grow- 
ing cluster ensemble than in the static one. This tends 
to reduce both the fractal dimension df and the walk 
dimension d^. If the effect is more pronounced for d^, 
then it may explain why an apparent measurement of 
the spectral dimension ds = 2df/dw was larger for FC in 
the growing cluster ensemble. Since the search algorithm 
is repeated to obtain an IFC cluster, the effect will be 
even larger in IFC than in FC percolation, as well as it 
is larger in FC than in ordinary percolation. 



V. SUMMARY 

In summary, we have introduced and studied an ex- 
tension to ordinary site percolation on a square lattice 
which we call iterated fully coordinated percolation. The 
static and dynamic critical behaviors of this model, as 
well as that of the ordinary and fully coordinated perco- 
lation models, were determined using methods such as 
Monte Carlo simulations, finite-size scaling, and normal 
mode analyses. All clusters used in our simulations were 
generated statically in grids of predetermined sizes on a 
square lattice, as opposed to growing them dynamically 
from a seed according to a dynamic rule. We found that 



the sizes of subregions of full connectivities in an incipi- 
ent infinite cluster greatly increases as we simulate from 
ordinary to fully coordinated and then to iterated fully 
coordinated percolation. The length scales of such subre- 
gions increase by a factor of two to three with each iter- 
ation. Although this implies a large increase in the local 
connectivities, our results also indicate that the ordinary, 
fully coordinated, and iterated fully coordinated perco- 
lation models have identical static and dynamic critical 
exponents within our statistical errors and thus we be- 
lieve that they belong to the same universality class, both 
statically and dynamically. 
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